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<N : Abstract 



^ I In this paper a detailed standard model (SM) calculation of the 0{as) virtual 

^ ' corrections to the decay width T{b —>■ sg) is presented {g denotes a gluon). 

Also the complete expressions for the corresponding 0{as) bremsstrahlung 

^Q . corrections to b ^^ sg are given. The combined result is free of infrared and 

^^ i collinear singularities, in accordance with the KLN theorem. Taking into 

^p I account the existing next-to-leading logarithmic (NLL) result for the Wilson 

r^ ' coefficient Cf , a complete NLL result for the branching ratio B {b — > sg) is 

II ■ derived. Numerically, we obtain B = (5.0ibl.O) x 10~^, which is more than 

^H. a factor of two larger than the leading logarithmic result B^^ = (2.2 it 0.8) x 

r~| , 10"'^. The impact of these corrections on the inclusive charmless hadronic 

branching ratio B^ of i?-mesons, which can be used to extract iKib/^bl in the 

k^ ■ context of the SM, is shown to be of similar importance as NLL corrections to 

;_( ■ 6-quark decay modes with three quarks in the final state. Finally, the impact 

of the NLL corrections to b ^ sg on B^ is investigated in scenarios, where 

the Wilson coefficient Cg is enhanced by new physics. 
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I. INTRODUCTION 

The theoretical predictions for inclusive decay rates of 5-niesons rest on solid grounds 
due to the fact that these rates can be systematically expanded in powers of AqcD/^^b [Sil' 
where the leading term corresponds to the decay width of the underlying 6-quark decay. As 
the power corrections only start at ©(Aq^jj^/m^), they affect these rate by at most a few 
percent. Theoretically, spectator effects of order 167r^(AQCD/'^6)^ [§,3 could be larger [Q, 
but for the decay rates of B^ and 5° they are experimentally known to be at the percent 
level as well [|]. Thus the accuracy of the theoretical predictions is mainly controlled by our 
knowledge of the perturbative corrections to the free quark decays. 

The inclusive charmless hadronic decays, B — >■ X^, where X^ denotes any hadronic 
charmless final state, are an interesting subclass of the decays mentioned above; as pointed 
out in ref. [§], a measurement of the corresponding branching ratio would allow the extraction 
of the presently poorly known ratio |Vufe/Vcf,|, where Vub and Vcb are elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. At the quark level, there are decay modes with three- 
body final states, viz. b —>■ q'q'q, {q' = u,d,s; q = d, s) and the modes b —>■ qg, with two-body 
final state topology, which contribute to the charmless decay width at leading logarithmic 
(LL) accuracy. Calculations of next-to-leading logarithmic (NLL) corrections to the three- 
body decay modes were started already some time ago in ref. , where radiative corrections 
to the current-current diagrams of the operators Oi and O2 were calculated, together with 
NLL corrections to the Wilson coefficients. Later, Lenz. et al. included in a first step |§] 
the contributions of the penguin diagrams associated with the operators Oi and O2, and in 
a second step ^ the same authors also included one-loop penguin diagrams of the penguin 
operators 03,...,0q; also the effects of the chromomagnetic operator Og were taken into 
account to the relevant precision needed for a NLL calculation. Up to contributions from 
current-current type corrections to the penguin operators, the NLL calculation for the three 
quark final states is now complete. 

In the numerical evaluations of the charmless hadronic branching ratio, the two body 
decay modes b —>■ qg were added in refs. [^Q at the LL precision, as the full NLL predictions 
were missing. To fill this gap, we recently wrote a short letter where NLL results for the 
branching ratio B{b -^ sg) were presented Q]. In the present work, we describe in detail the 
non-trivial two-loop NLL calculation, which led to the results in ||^. As the NLL corrections 
enhance B{b -^ sg) by more than a factor of 2, we also analyze in the present paper their 
impact on the charmless hadronic branching ratio. 

The decay b ^ sg gained a lot of attention in the last years. For a long time the 
theoretical predictions for both, the inclusive semileptonic branching ratio Es\ and the charm 



multiplicity n^ in S-meson decays were considerably higher than the experimental values [|10 
An attractive hypothesis, which would move the theoretical predictions for both observables 
into the direction favored by the experiments, assumed the rare decay mode b ^ sg to be 
enhanced by new physics. 

After the inclusion of the complete NLL corrections to the decay modes b — >■ cuq and 
b -^ acq {q = d,s) [jll|, the theoretical prediction for the semileptonic branching ratio and 



the charm multiplicity [§] are 



;B*I^ = (11.7 ±1.4 ±1.0)% , n**^ = 1.20 ± 0.06 , (1) 

where the second error in Blf takes into account the spectator effects estimated in ref. 0. 
The experimental results from measurements at the T(4S) resonance and those from the Z^ 
resonance at LEP and SLD were recently summarized [|12| to be 

i3j(^^^ = (10.45 ± 0.21)% , nj^^"^^ = 1.14 ± 0.06 , 

Sf = (10.79 ± 0.25)% , nf = 1.17 ± 0.04 . (2) 

We would like to stress that in the theoretical results the renormalization scale was taken in 
the interval [mb/4, 2m;,]. If one only considers /i G [mf,/2,2mb], the theoretical predictions 
would only have marginal overlap with experimental data. This implies that there is still 
room for enhanced b ^ sg. We therefore also illustrate in this paper the influence of the 
NLL corrections to b -^ sg on the charmless hadronic branching ratio in scenarios where the 
Wilson coefficient Cg is enhanced by new physics. 

We also would like to mention that the component b -^ sg oi the charmless hadronic 
decays is expected to manifest itself in kaons with high momenta (of order mb/2), due to its 
two body nature [0]. Some indications for enhanced b —>■ sg in this context were reported 



by the SLD collaboration |jT^. For a review of other hints for enhanced b —>■ sg, the reader 
is referred to [p^ . 

Within the SM, the LL prediction for the branching for b ^ sg is known to be B{b — *■ 
sg) ~ 0.2% H^. The process b —>■ sgg, which gives a NLL contribution to the inclusive 



charmless decay width has already been studied in the hterature ||T7| . |T8|| . In [T^ a complete 
calculation was performed in regions of the phase space which are free of collinear an infrared 
singularities. Putting suitable cuts, the branching ratio for b -^ sgg was found to be of the 
order 10~^ in these phase space regions. A complete calculation requires the calculation 
of a regularized version for the decay width r(6 -^ sgg) in which infrared and collinear 
singularities become manifest. Only after adding the virtually corrected decay width r(6 — >■ 
sg) a meaningful physical result is obtained. In addition, as we will see later, also the tree 
level contribution of the operator Og to the decays b —>■ sff, with f = u,d, s, has to be 
included. 

The remainder of this paper is organized as follows: In section |I|, we review the theoret- 
ical framework and discuss the steps needed for a NLL calculation for B{b —>■ sg). Section 
pn| is devoted to the calculation of the virtual corrections to the matrix elements {sg\0i^2\b), 
including renormalization, while section |IV| deals with virtual corrections to {sg\Os\b). In 
section ^ the virtual corrections to the decay width T(b ^ sg) are calculated. In sec- 
tion |VI| the decay width T(b ^ sff) is given. Sections |VII| and |VII1| deal with the gluon 
bremsstrahlung matrix elements {sgg\0i^2.s\b) and the corresponding decay width, respec- 
tively. The analytic results for the NLL branching ratio B{b — > sg) can be found in section 
p(| , while the numerical evaluations are presented in section ^. Section ^ deals with the 



impact of the NLL corrections to B{b — * sg) on the charmless hadronic branching ratio in 
the standard model, while in section PCI1| similar questions are addressed in scenarios where 
the Wilson coefficient Cg is enhanced by new physics. We conclude with a short summary 
in section PCII1| . An explicit parametrization of the NLL Wilson coefficient C^{mb) is given 
in appendix |^. 



II. THE EFFECTIVE HAMILTONIAN 

We use the framework of an effective low-energy theory with five quarks, obtained by in- 
tegrating out the heavy degrees of freedom, which in the SM are the t-quark and the 
PF-boson. We take into account operators up to dimension six and we put m^ = 0. In this 
approximation the effective Hamiltonian relevant for radiative decays and b -^ sg{g) reads 

AC ^ 

Kff = ^ V,lVt, Y: C^if^Wfi), (3) 

where Gp is the Fermi coupling constant and Ci{^) are the Wilson coefficients evaluated 
at the scale /i; Vtb and Vts are matrix elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. The operators Oj read [jl9| 



O5 = {sLlf^lulpbL) E.iqrrYq) , O, = {slI.I.IpTHl) E^iq^YYT^'q) , 



flV 



In the dipole operators 0-j (Os), e and F^jy {qs and G"^^) denote the electromagnetic (strong) 
coupling constant and field strength tensor, respectively. T"^ [A = 1, ...,8) are SU{3) color 
generators; L = (1 — 75)/2 and R = {1 + 75)/2 stand for left- and right-handed projectors. 
In eq. (^), mb(/i) is the running 6-quark mass in the MS scheme at the renormalization scale 
fi. Henceforth, mq{fi) and rriq denote MS running and pole masses, respectively. To first 
order in a^, these masses are related through: 

It is well-known that QCD corrections to the decay rate for 6 — ;> 57 bring in logarithms 
of the mass ratios rrih/rnw and rrib/mt. The same is true for the process b ^ sg: QCD 
corrections to this process induce terms of the form asa"ln™'(?7ib/M), where M = rrit or 
my/ and m < n (with m,n = 0,1,2, .. .). 

One can systematically resum these large terms by renormalization group techniques. 
Usually, one matches the full standard model theory with the effective theory at a scale of 
order ruw At this scale, the large logarithms generated by matrix elements in the effective 
theory are the same ones as in the full theory. Consequently, the Wilson coefficients only pick 
up formally small QCD corrections. Using the renormalization group equation, the Wilson 
coefficients are then calculated at the scale /x = /i;, ~ m^, at which the large logarithms are 
contained in the Wilson coefficients, while the matrix elements of the operators are free of 
them. 

So far the decay rate for b ^ sg has been systematically calculated only to leading 
logarithmic (LL) accuracy, i.e., for m = n. 

A consistent calculation for b —>■ sg at LL precision requires the following steps: 



1) the extraction of the Wilson coefficients from a matching calculation of the full stan- 
dard model theory with the effective theory at the scale /i = fxw to order a°; fi]y 
denotes a scale of order mw or m^; 

2) a renormalization group treatment of the Wilson coefficients, using the anomalous- 
dimension matrix to order a]; 

3) the calculation of the decay matrix elements {sg\CiOi\b) at the scale n = jii,to order (yf^; 
/ife denotes a scale of order m},. We note that the matrix elements associated with the 
four Fermi operators (z = 1 — 6) can be absorbed into the effective Wilson coefficient 
C|^, when working at LL precision. In the naive dimensional regularization scheme 



(NDR), which we use in this paper, one obtains ||T9| 



Cf = Q + C3-|C4 + 20C5-f Ce . (6) 

From the analogous decay 6 — * 57 it is well-known that next-to-leading logarithmic 
(NLL) corrections drastically reduce the large renormalization scale dependence of the LL 
branching ratio. This implies, in particular, that the NLL corrections are relatively large, at 
least for certain scales (within the usually considered range mb/2 < fit < ^nib). Motivated 
by the situation in this analogous process, we present in this paper a systematic calculation 
of the NLL corrections to b ^ sg. 

The principal organization of such a calculation is straightforward: Each of the three 
steps listed above has to be improved by going to the next order in a^: 1) The matching 
has to be calculated including ag corrections; 2) the renormalization group treatment of the 
Wilson coefficients has to be performed using the anomalous dimension matrix to order a^; 
3) finally, the order as corrections to the decay matrix elements have to be worked out. We 
note that this step involves both, the calculation of virtual- and hremsstrahlung corrections 
to 6 — >• sg. 

The ffist two steps are already available in the literature. The order a^ matching of 
the dipole operators O7 and O^ was calculated in refs. |2^, while the matching conditions 
and the anomalous dimension matrix for the four Fermi operators have been calculated by 
several groups [^]. These calculations were done in the "old operator basis", introduced by 
Grinstein et al. |2^. The most difficult part, the order a^ mixing of the four- Fermi operators 



into the dipole operators requires the calculation of three loop diagrams [|l9l. In order to 
perform a consistent NDR calculation (i.e. with anticommuting 75), the old operator basis 
was replaced by the new one displayed in eq. @) . The full 8x8 anomalous dimension matrix, 
the corresponding matching conditions and the definition of the evanescent operators is given 
in ref. [1^ and is repeated in appendix ^ of the present paper. 



Step 3), the calculation of the virtual 0{as) corrections to the matrix elements Mj = 
{sg\Oi\b), as well as the evaluation of the gluon bremsstrahlung process h -^ sgg, is performed 
the ffist time in the present paper. As illustrated in table |, the LL Wilson coefficients 
C^{fib),...,CQ{fib) are much smaller than C^^fib) and C2{fib)- We therefore only calculate 
Ml, M2 and Mg together with the corresponding bremsstrahlung corrections. As Mi and 
M2 vanish at one-loop (i.e. without QCD corrections), only the leading order pieces, C°(y^b) 
and 6*2 (/Xfe), appearing in the decomposition 



Qifi,) = C°(/i,) + 



Air 



Clif^t 



(7) 



of the NLL Wilson coefficients Ci{fib) and C2(/Uf,) are needed. On the other hand, the 
operator Og contributes to Mg already at tree-level. Consequently the full NLL Wilson 
coefficient C^^{fiiy) is needed. The numerical value of the NLL piece Cg'^ (defined as in eq. 
(^) and (|^)) is also given in table |, while the analytic form is relegated to appendix ^ 





jJL = mw 


fi = 9.6 GeV 


/x = 4.8 GeV 


H = 2.4 GeV 


as 


0.121 


0.182 


0.218 


0.271 


C? 


0.0 


-0.335 


-0.497 


-0.711 


C° 


1.0 


1.012 


1.025 


1.048 


CI 


0.0 


-0.002 


-0.005 


-0.010 


cl 


0.0 


-0.042 


-0.067 


-0.103 


c? 


0.0 


0.0002 


0.0005 


0.001 


Q° 


0.0 


0.0005 


0.001 


0.002 


CO 


-0.192 


-0.285 


-0.324 


-0.371 


cl 


-0.096 


-0.136 


-0.150 


-0.166 


^0,eii 


-0.196 


-0.280 


-0.314 


-0.356 


<^8 


-0.097 


-0.135 


-0.149 


-0.165 




-2.166 


-1.318 


-1.098 


-0.950 


r'eff 
<^8 


-0.118 


-0.154 


-0.168 


-0.186 



TABLE I. Wilson coefficients CP(^) {i = 1,...,8), Cl'''^ and Cf (see eq. (0) in the text) at 
the matching scale ^i = raw = 80.33 GeV and at three other scales, ^ = 9.6 GeV, ^ = 4.8 GeV 
and /i = 2.4 GeV. For as(//) (in the MS scheme) we used the two- loop expression with 5 flavors 
and as{mz) = 0.119. The entries correspond to the pole top quark mass mt = 175 GeV. 



III. VIRTUAL CORRECTIONS TO Oi AND O2 



In this section we present the calculation of the matrix elements of the operator Oi and 
O2 for h ^ sg xvp to order a^ in the NDR scheme. The one-loop (a^) matrix elements 
vanish and we must consider several two-loop contributions. Since they involve ultraviolet 
singularities also counterterm contributions are needed. These are easy to obtain, because 
the operator renormalization constants Zij are known with enough accuracy from the order 
as anomalous dimension matrix |19(| . 



A. Regularized two-loop matrix elements of 0\ and O2 



For the following discussion it is useful to define the operators Oi and O2: 




r'=* ~rTF ~rTT 




0" °V,^t^ °V_>6%, '^'"'"- 

a) c) e) g) 

b s b s b s b 





%^D' %LJ'' J'^^^^"" '^^ 

b) d) ^ Q h) 

FIG. 1. Graphs associated with the operators Oi and O2. The wavy hnes represent gluons; the 
real gluons are understood to be attached to the circle-crosses. 
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FIG. 2. Building block /^ (with an off-shell gluon) for the diagrams a), b), e) and f) in fig. |l| 
and building block J^fs for the diagrams c) and d) in fig. 1^. g* and g denote an off-shell and an 
on-shell gluon, respectively. 

Oi = 20i + i02 ; 02 = 02 . (8) 



Oi and O2 are nothing but the current- current operators in the old basis p^ : 

Ol = (sLalt^CLp) (cLpYha) , O2 = (sLp'J^.CLp) {cial^^La) ■ (9) 

We now present the calculation of the matrix elements Mj = {sg\Oi\h)\ The dimensionally 
regularized matrix element M2 is obtained by calculating the two-loop diagrams a) - h) 
shown in fig. |I]. 

We start with the calculation of the diagrams a) - f ) in fig. |l|, in which the virtual 
gluon connects the charm quark in the loop with an external fermion legQ. The main steps 
of the calculation are the following: We first calculate the Fermion loops in the individual 
diagrams, i.e., the 'building blocks' Ip and J^p shown in fig. 0; J^p denotes the sum of the 
the two diagrams on the right. 

We work in d = 4 — 2e dimensions; the results of the building blocks are presented 
as integrals over Feynman parameters after integrating over the (shifted) loop-momentum. 



^ The diagrams g) and h) are much easier to calculate than those in a) ~ f), because nic is the 
only scale in the corresponding integrals. 

6 



Then we insert these building blocks into the full two-loop diagrams. Using one more 
Feynman parametrization, we calculate the integral over the second loop-momentum. As 
the remaining Feynman parameter integrals contain rather complicated denominators, we 
do not evaluate them directly. At this level we also do not expand in the regulator e. The 
heart of our procedure which will be explained more explicitly below, is to represent these 
denominators as complex Mellin-Barnes integrals p3|. After inserting this representation 
and interchanging the order of integration, the Feynman parameter integrals are reduced to 
well-known Euler Beta-functions. Finally, the residue theorem allows to write the result of 
the remaining complex integral as the sum over the residues taken at the pole positions of 
Beta- and Gamma-functions; this naturally leads to an expansion in the ratio z = {mc/mb)'^, 
which numerically is about z = 0.1. 

We express the diagram on the left in fig. ^ (denoted by /^) in a way convenient for 
inserting into the two-loop diagrams. As we will use MS subtraction later on, we introduce 
the renormalization scale in the form /i^exp(7£;)/(47r), where 7^ ~ 0.577 is the Euler con- 
stant. Then, MS corresponds to subtracting the poles in e. In the NDR scheme, J/3 is given 
by2 



9s 



r(e) 12^' exp(7£;e) (1 - e) exp(i7re) T^ (r/3/ - r^-ff^) L x 



47r2 



ni'i 



x{l — x) 



+ i5 



(10) 



where r is the four-momentum of the (off-shell) gluon, rric is the mass of the charm quark 
propagating in the loop and the term i5 is the " e-prescription" . The free index (3 will be 
contracted with the gluon propagator when inserting the building block into the two-loop 
diagrams a), b), e) and f) in fig. |1]. Note that Ip is gauge invariant in the sense that r^Ip = 0. 

Next we give the sum of the two diagrams on the right in fig. 0, using the decomposition 
^ . The on-shell gluon has momentum q, color A and polarization a (therefore we drop 



m 



the terms q and qa), while the off-shell gluon has momentum r, color B and polarization (3. 
This building block, denoted by J^/f 5 can be decomposed with respect to the color structure 
as 



jAB 

■^0/3 



T+(g,r){T^r^}+T-^(g,r) 



rpA^jt 



(11) 



The quantities T\{q,r) and TJq.r) read 



ap\ 



TaaiQ^r) 



327r2 



E{a, P, r) Ais + E{a, p, q) Aig - E{I3, r, q) -^Aiss 

(gr) 

-E{a, r, q) T-^^^25 - E{a, r, q) T^ ^^26 



(12) 



^The fermion/gluon and the fermion/photon couplings are defined according to the covariant 
derivative D = d + igsT^A^ + ieQA where T^ 



— are the SU(3) generators. 



Taail^r) 



327r2 



gaf3 Az2 + ^ gap Ais + 7^ r„ Aig + 7a ^/3 Azii + 7^, g^g Azi2 

+/ — — A2i5 + f — — A2i7 + ^ — — A2i9 + ^ — — A«21 1^ • 

(grj (gr) [qr) [qr) 



The matrix £" in eq. (jT^) is defined as 

E{a, /3, r) = 7^,7/3/ - 7ar/3 + 7/3(r„) 



Jal3- 



(13) 



(14) 



In a four- dimensional context these E quantities can be reduced to expressions involving the 
Levi-Civita tensor, i.e., E{a,l3,'~f) = —"^ £^0/37^ 7^75 (in the Bjorken-Drell convention). The 
dimensionally regularized expressions for the Az functions read 



Ai2 = AB dxdyC ^ "il- x) \Aiqr)xyt -2r'^x^t + r^xe + C 

—2r'^x'^ye + 2r^x^e + r^xye — 2r'^xe + yC — C] 
Ais =—45"/ dx dy C~^^''[A{qr)x^ye + 2{qr)xye 

-2r'^x^e + r^x^e + r^xe + xC + C] 



-f-l-e/ 



Aiu= AB- dxdyC-\l-x) 



A(qr)xye — 2{qr)xe — 2rxe + rxe + C 



Aii2 = AB / dxdyC ^ ''[A{qr)xy'^e + 2{qr)xye 

—2r'^x^yt — 2r^x^e + r^xyt + 2r^xe + yC + C\ 
Aii5 = l6B-{qr)e f dxdyC~^-'x^{l-x) 

Aii7 = -8B-{qr)e f dxdyC'^''xy{l-2x) 

Aiig = 8B^{qr)e dxdyC^^"^x{\ — x — y + 2xy) 

Ai2i = 8B-(qr)e f dxdyC-^-'xyil - 2y) 
Js 



(15) 



Ais = —AB^ / dxdyC ^ ^[A{qr)x'^yt — A{qr)xyt 

-2r^x^e + Sr^x^e - r^xe + 3xC - C] 

Aig = AB^ dxdyC~^~^[A(qr)xy'^e — A(qr)xye 
Js 

-2r'^x^ye + 2r^x^e + r^xye - 2r^xe + 3yC - C] (16) 

Ai23 = -Ai26 = 8S+(gr)e f dxdyC~^-'xy (17) 

Ai25 = -8 5+(gr)e f dx dyC-^-'x{l - x) (18) 



(19) 

(20) 

(21) 
(22) 

(23) 
(24) 

(25) 

(26) 

(27) 



where C, C and B are given by 



C = m^ — 2xy{qr) — x(l — x)r 
C-^^' = - exp(me) [a;(l - x)]-^-' 



'-id 






I- X 


2 

x[l — x) 


-1-e 



(28) 



B+ = (l + e)r(e)exp(7se)/i2% B" = (e - l)r(e) exp(7se)/i2^ . (29) 

The range of integration in (x, y) is restricted to the simplex S", i.e., < y < (1 — x) and 
< X < 1. 

We are now ready to evaluate the two-loop diagrams. Due to the absence of extra 
singularities in the limit of vanishing strange quark mass, we set m^ = from the very 
beginning. 

In ref. [Q the detailed calculation of one of the diagrams in fig. |l|a) was presented for 



h -^ S7. As all the other diagrams, which involve the building block Ip, i.e., a), b), e) 
and f) in fig. |1|, can be computed in a very similar way, we prefer to concentrate on the 
diagrams involving the building block J^p. As an example in this class, we concentrate on 
the diagrams d) in fig. |1|, which we redisplay in fig. |^ in order to set up the notation for the 
momenta. 



b(p) p+r s(p') 



"^^ 







FIG. 3. Feynman diagram for the Mellin-Barnes example. The momentum and the polarization 
vector of the emitted gluon are denoted by q and e, respectively. 

The sum M2{d) of the two diagrams can be decomposed into a color symmetric part 
M2 (d) and a color antisymmetric part Mg" (d) according to 

M^(d) = M+{d) + M^{d) , (30) 

with 

M,(.)^..(-.)/-T-T-.^^^i/|^.(p') {T,,s^) ^^/^(.)^ 



Mm - 0. Y "■^'iS 7 / l^^f-'' (^- ^") 7?^ -* "(rt^ ■ (31) 



where T^^ and T„^ are given in eqs. ([121) and ([131) , respectively. As the calculation of M^id) 
is nothing but a repetition of the b ^ sj case, we concentrate on M2 (d) in the following. 
All the Ai quantities in T~g contain the factor C"^""^, whose explicit form is given in eq. 



(28). M2 (d) can be written in the form 



iJs ( ■\ fABCrpBrpC 



M-(d) = ^^ (-i)f^^^T^T 



327r2 



with Di = (r^ + 2(pr)), D2 = r^, D3 = r'^ + 2{qr)y/{l — x) — ml/{x{l — x)). The symbol P(r) 
is a matrix in Dirac space, which depends in a poljTiomial way on the integration variable 
r. In the next step, the three propagators Di, D2 and D^ in the denominator are Feynman 
parametrized as 



1 r(3 + e) f dudww'' 



D1D2DI+' r(l + e) is [r2 + 2{pr)u + 2{qr)yw/{l - x) - m>/(x(l - x)) + tSf'^' 



(33) 



with < w < 1 — u and < m < 1. Then the integral over the loop momentum r is 
performed. At this level, a four dimensional integral over the Feynman parameters {x, y; u, w) 
remains. It is useful for the following to perform the substitutions 

/ (1 -a;')(l -W -y') n '\ ' ' ' f'jA\ 

x^x] y^ ; u^[\ — w)u] w^uw. (34) 

The new variables then run in the intervals 

x,m,m;'g [0,1]; ye[l-w\l\. (35) 

Taking into account the corresponding Jacobian and omitting the primes (') of the integra- 
tion variables, M2" {d) can be cast into the form 



M2{d) = ^,{-i)f^^''T^T<' 



dxdydudw u{p') 



' C ^ 1 ^ 1 



(J2e (72e '^ (^l+2e 



uip) , (36) 



where Fi, F2 and F3 are matrices in Dirac space depending on the Feynman parameters x, 
y, u, w. Note that this expression is understood to be written in such a way that Fi, F2 and 
F3 are independent of rric. The charm quark mass then only enters through C, which reads 
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C = miuy{l-w)+ ' w . (37) 

X{1 — X) 

In what follows, the ultraviolet e regulator remains a fixed, small positive number. 

The central point of our procedure is to use now the Mellin-Barnes representation of the 
denominators that look like propagators {l/{k^ — M^)'^) ||23|, which is given by (A > 0) 



(p - M2)^ (p)^ r(A) 

10 



The symbol 7 denotes the integration path which is parallel to the imaginary axis (in the 
complex s-plane) hitting the real axis somewhere between —A and 0. In this formula, the 
"momentum squared" A;^ is understood to have a small positive imaginary part. 



In our approach, we use formula ( |5^ ) in order to simplify the remaining Feynman param- 
eter integrals in eq. ( pG] ) where we represent the factors l/C^^ and 1/C^^'^'^ as Mellin-Barnes 
integrals using the identifications 

k"^ ^ niluy{l - w) ■ M^^^^^ . (39) 

By interchanging the order of integration, we first carry out the integrals over the Feynman 
parameters for any given fixed value of s on the integration path 7. These integrals are 
basically the same as for the massless case rric = (in eqs. (|36| ) and (0)) up to the factor 



w 



uy [1 — w) x{l — x) 




(40) 



Note that the functions Fi, F2 and F3 are such that the Feynman parameter integrals exist 
if the integration path 7 is properly chosen. In the terms involving F2 and F3 in eq. (p6|), 
the path must be chosen such that — e < Re(s) < 0; in the terms involving Fi the situation 
is slightly more complicated: C in the numerator should be replaced by the r.h.s of eq. (p7|). 
For the terms proportional to ml the path has to be chosen as for the F2 and F3 contributions. 
The terms proportional to m^, however, lead to Feynman parameter integrals which do not 
converge for values of s on this path. It turns out that the path has to be chosen such that 
— 2e < Re(s) < — e in order to have convergent integrals for these terms. 

We would like to mention that the variable substitutions in eq. ( |3^ were constructed in 
such a way that all the Feynman parameter integrals are either elementary or of the form 

/o dx xP{l - xY = I3{p + l,q + I). 

For the s integration we use the residue theorem after closing the integration path in the 
right s-halfplane. According to the above discussion, the residue at s = — e has to be taken 
into account in the terms proportional to ml. In the other terms, however, the residue at 
s = — e must not be taken into account. The other poles inside the integration contour are 
located at 

5 = 0,1,2,3, 

s = 1 — e, 2 — e, 3 — e, 

s = 1 -2e,2-2e,3-2e, 

s = 1/2 - 2e, 3/2 - 2e, 5/2 - 2e, 

s = l-3e,2-3e,3-3e, (41) 

The other two-loop diagrams are evaluated similarly. The non-trivial Feynman integrals 
can always be reduced to /5-functions after suitable substitutions. 

The sum over the residues naturally leads to an expansion in z = {ml /ml) through the 
factor [m^^/m^lY in eq. (|3)- This expansion, however, is not a Taylor series; it also involves 
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logarithms of z, which are generated by the expansion in e. A generic diagram which we 
denote by G has then the form 



7^ 

G' = co + Ecnn^^"ln'";^ , z = ^ . (42) 



n.m 



The power n in eq. (^) is in general a natural multiple of 1/2 and m is a natural number 
including 0. In the explicit calculation, the lowest n turns out to be n = 1. This implies the 
important fact that the limit rric -^ exists. 

From the structure of the poles one can see that the power m of the logarithm is bounded 
by 4, independent of the value of n. For a detailed explanation, we refer to [^]. As in this 
reference, we retain all terms up to ra = 3 in our results. 

Unlike in 6 —» 57, the diagrams in the individual figures are not gauge invariant. This 
statement holds even for the sum of all the diagrams in a) - f ) in fig. |I|. A gauge invariant 
result is only obtained after including the diagrams in g) and h)[|. We would like to mention 
that the diagrams analogous to g) also exist for b — > 57. Their sum, however, vanishes 
in this case. As there are no gauge invariant subsets, we only present the result which is 
obtained by summing all diagrams a) - h) in fig. |I|. The result for M2 = {sg\02\b) reads 
(using z = {nic/mbY and L = In^): 

2592 71 \ f^ J 

384 

{ 2170 - 547r2 + ^[48816 - 25271^ + (22680 - 16207i'^)L 

t 

+2808L2 + 612L3 - 6480C(3)] 

-1267271^^2/2 + z^ [66339 + 18727r2 + (-40446 + 15127r2)L 

+6642L2 - 1008^2 + 7776C(3)] 

+^^[-3420 - 607r2 - 6456L + ISML'^] 

+2Am[ - 28 + z(549 - 247r2 + 153L + 72L^) 

+^^(-432 + 307r2 + 54L - 90L2) + z^{-2hQ + 192L)]} . (43) 

In this expression, the symbol C, denotes the Riemann Zeta function, with ^(3) ~ 1.2021; 
The symbol {sg\0^\b)tj.ce denotes the tree level matrix element of the operator Og. As such, it 
contains the running 6-quark mass and the running strong coupling constant, both evaluated 
at the scale // (see eq. @)). However, as the corrections to O2 are explicitly proportional 
to Usi we are allowed (modulo higher order terms) to identify the running 6-quark mass 
with the pole mass mb'-, in the same spirit we can identify the strong coupling constant with 
gsijni,). With this interpretation, which we will use in the following, {sg\0^\b)t^ee is a scale 
independent quantity, reading 



^We thank M. Neubert for making us aware of these diagrams. 
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[sg\Os\b)tree = m,^^u{p')^<iRT^u{p) . (44) 



We now turn to the matrix elements of the operator Oi. Due to the specific color 
structure it is straightforward to see that only the diagrams e) and f) in fig. |l| yield a 
non- vanishing contribution, which is generated by the color symmetric part of the building 
block Jap in eq. ([TT|) . The complete regularized result for Mi = {sg\Oi\b) reads 

1 « / X -4e 

^1 = 7^ — {s9\Os\b)tree 

yo IT 

{ - — - 87 + ^[120 - levr^ + (120 - 367r2)L 

+12L^ + 4L^ - 144C(3)] 

+z'^[84: + 327r^ -247T^L 

-12L'^ + iL^] + z^[-56 - 127r2 + 96L - 36L^] 

-4:ni[3 + ^(-24 + 271^ - 6L - 6L^) 

+z'^{-6 + 2n^ + l2L~6L'^)-l2z^]}. (45) 

The regularized matrix elements Mi and M2 of Oi and O2 in the operator basis (^) are 
related to Mi in eq. (|45|) and M2 in eq. (|43|) as follows: 



Ml = i Ml - I M2 ; M2 = M2 . (46) 



B. Counterterms to the Oi and O2 contributions 

The operators mix under renormalization and thus the counterterm contributions must 
be taken into account. As we are interested in this section in contributions to b ^ sg 
which are proportional to Ci and C2, we have to include, in addition to the two- loop matrix 
elements of CiOi and C2O2, also the one-loop matrix elements of the four Fermi operators 
CiSZijOj {i = 1,2; j = 1,...,6) and the tree level contribution of the magnetic operator 
CiSZisOs {i = 1,2). In the NDR scheme the only non- vanishing contributions to b —^ sg 
come from j = 4, 8 only. The operator renormalization constants Zij are obtained from the 
leading order anomalous dimension matrix in the literature |19[ Q. The entries needed in our 
calculation are 



^Note that the effective anomalous dimension matrix ^*''™ given in [19| has to be converted into 
7*^, before the relevant (5Z-factors can be read off. 
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SZ,, = ^ , 5Z,s = ^ . (48) 

Dvre 1087re 

The counterterm contributions Mf' and M|* proportional to Ci and C2 are then given by 
Mf = (.sI^ZmO, + SZMb) = (^ 1 (!;;£)"■ + ^ J^ 1] (.,,|0«|6).„ . (49) 



Mf = (891^^2404 + SZ^Oilb) = I -^ i ('^] +-^-\ {sg\0,\b)„„ . (60) 

doTT e \ n I TT 108 e 



We note that there are no one-loop contributions to the matrix element for b -^ sg from 
the counterterms proportional to the evanescent operators Pn and P12 given in appendix A 



of ref. 19 . 



C. Renormalized matrix elements of Oi and O2 

Adding the regularized two- loop result in eq. ( ^3D and the counterterm in eq. (|5DD, we 
find the renormalized result for M2 in the NDR scheme: 

Mr = {sg\Os\b)tree ^ U In ^ + v^) , (51) 



with 



(52) 



^ 27 

Re(r2) = { - 2170 - 5471^ + z[48816 - 25271^ + (22680 - 16207r^)L 

+2808L2 + 612^3 - 6480C(3)] 

-1267271^2^/2 + 22[66339 + 18727r2 + (-40446 + 15127r2)L 

-1-6642^2 - 1008L^ + 7776C(3)] 

+z^[-3A20 - 607r2 - 6456L + 7884^^]} 



Im(r2) = — { - 28 + ^[549 - 2471^ + 153L + 721"^] 

+z'^ [-432 + 307r2 + 54L - 90/.^] + z^ [-259 + 192L] } (53) 

Here, Re(r2) and Im(r2) denote the real and the imaginary part of r2, respectively. The 
quantity z is defined as z = {ml/ ml) and L = \n{z). 
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0.25 




9 / 9 

z = m^imi 

FIG. 4. Real and imaginary part of r2 in the NDR scheme (from eq. (p^)). 

Similarly, we obtain the renormalized version of Mi by adding the regularized two-loop 
result in eq. (H) and the counterterm in eq. (|^); we find 



Mr = {sg\0,\h),,,,^ Uln'^ + n 



(54) 



with 



173 
162 



(55) 



Re(ri) 



1 



{4877 - 547r2 + 36z[1086 + 297r2 + (360 + 367r2)L 
3888 

+51^2 + 8L3 + 144C(3)] 

-126727r2;z^/2 + '^z^[mib - 807?^ + (-4494 + 3847r2)L 

+864L2 - 148L3 + 864C(3)] 

+122^ [93 + 767r2 - 1186L + 900L2]} 



Ini(ri) 



TT 



324^ ^ ^ 

+Qz^[-in + 197r2 + 72L - 57L^] + 2z^[-421 + 192L]} 



(56) 



In figs. ^ and |] we show the real and the imaginary parts of r2 and ri, respectively. For 
2 > 1/4 the imaginary parts must vanish exactly; indeed we see from these plots that the 
imaginary parts based on the expansion retaining terms up to z^ indeed vanish at 2; = 1/4 
to high accuracy. 



15 



0.05 



0.10 



0.15 



0.20 



0.25 




z = ml I ml 

FIG. 5. Real and imaginary part of ri in the NDR scheme (from eq. (pq)). 

IV. VIRTUAL CORRECTIONS TO Os 

In this section we calculate the order a^ virtual corrections to the matrix element 

Ms = {sg\0,\h) . (57) 

As the contributing Feynman graphs in fig. |^ are one loop diagrams, the computation of 
Mg is straightforward. We use dimensional regularization for both, the ultraviolet and the 
infrared singularities. Singularities which appear in the situation where the virtual gluon 
becomes almost real and coUinear with the emitted gluon are also regulated dimensionally; 
on the other hand, those singularities where the almost real internal gluon is collinear with 
the s-quark, are regulated with a small strange quark mass tti^; the latter manifest themselves 
in logarithmic terms of the form ln(p), where p = {ms/mhY- 

We were able to separate the ultraviolet 1/e poles from those which are of infrared 
(and/or collinear) origin. For ultraviolet poles we use the symbol 1/e in the following, while 
collinear and infrared poles are denoted by 1/eiR. 

When working in Feynman gauge for the gluon propagator, the individual diagrams 
contributing to Mg have the following infrared/collinear properties (the letters refer to the 
diagrams in fig. ^: a) and b) are free of infrared and collinear singularities; c) has combined 
infrared/collinear singularities of the form 1/ejj^ or ln(p)/eiR as well as l/ejR poles; d) has 
combined infrared/collinear singularities of the form l/efj^ as well as 1/eiR poles; e) has a 
collinear singularity of the form ln(p); f) is free of infrared and collinear singularities; g) 
has a combined collinear and infrared singularity of the form ln(p)/eiR as well as collinear 
singularities of the form In^(p) and ln(p); h) has an infrared singularity of the form l/ei^; 
more precisely, this diagram is proportional to the combination (1/e — l/ei^). 
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FIG. 6. Diagrams associated with the operator Os- The real gluon can be attached to any of 
the circle-crosses on the fermion lines. 



As the results of the individual diagrams are not very instructive, we only give their sum: 



«., 



Ms = ^fs{sg\Os\b),,,,, 



(58) 



with 



/s 



3 _ (4 ln(p) + 9 + 9i7r) 11 



"IR 



3eiR 



1 
3 



597r2 
^2" 






-2e 



+ 



+ 1 - 8\n{p) + 2W{p) - 8in 



(59) 



We would like to mention that we did not include diagrams with self energy insertions in 
the external legs. As we work in an on-shell renormalization scheme with respect to quark 
and gluon fields, such diagrams are cancelled against counterterm contributions. 



A. Counterterms to the Og contribution 



The counterterm is generated by expressing the bare quantities in the tree-level matrix 
element of Os by their renormalized counterparts. It has the structure 



Mi' = 6R{sg\0s\b)tree , 



(60) 



where the factor 6R is given by 



6R = JZ2{mb) JZ2{ms) J Z^ Zg^ Z^, ^88 - 1 



(61) 



Zi2{f^b)i Z2{ms) and Z^ denote the on-shell wave function renormalization factors of the b- 
quark, the s-quark and the gluon, respectively. Zg^ and Z^^ denote the MS renormalization 
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constants for the strong coupling constant Qs and the 6-quark mass factor, which appear 
exphcitly in the definition of the operators (see eq. (D). Finally, Zgs is the renormalization 
factor of the operator Og. 



The explicit form of Z2{m) reads 






-2e 



1 2 

- + — + 4 
e eiR 



(62) 



where we again separated infrared and ultraviolet poles. For Z^ we get in the on-shell 
scheme: 



1 + 



as 5 /I 
2^2 



1 
eiR 



Os 1 
27r3 



E 



2 In 



/i 



(63) 



The sum in this formula run over the five flavors f = u,d, c, s, b. For Z^^ and Zgg (see ref. 
iTOll) one obtains 



^m,j 



Os 4 

Air e 



1 + 



Att 3e 



(64) 



Finally, the renormalization constant for the strong coupling constant reads 

a. 



1- 



Att 



2 3 



A^ 



/ 



(65) 



Inserting the various Z factors in eq. (|6T|) , one obtains 



An 



11 31 , rrih 2 .5-^, rrit 16 

— + - 81n— --Vln^^ + — -21np 

H 6 ^ n 3 



3e 6eiR 



(66) 



B. Renormalized matrix element of Og 



Adding the regularized matrix element of Og in eq. (|58D and the counterterm contribution 
Mg* in eq. (|60D, one obtains the renormalized result 



a.. 



Mr=^fr{s9\0s\b)tree 



(67) 



with 



/s 



3 (8 ln(p) + 49 + ISivr) 



-IR 



6eiR 






-2e 



29, mb 

— m — 

3 /i 



m/ 



+^Ein^-5 + 



/^ 



597r2 
~36" 



- glnp+gln p 



-trr 



(68) 



In eq. ( |68D the sum runs over the five fiavors / = u,d,c,s,b, and p = (mg/mb)^. We 
anticipate that the singular terms of the form l/ef^, l/ejR and Inp in eq. ( |68|) will cancel 
(at the level of the decay width) against the corresponding singularities present in the gluon 
bremsstrahlung corrections to 6 — >■ sg. On the other hand, the logarithmic terms ln{mf/p), 
which also represent some kind of singularities for the light fiavor f = u,d,s are not cancelled 
by the gluon bremsstrahlung process. Keeping in mind that these terms originate from the 
renormalization factor Z3 of the gluon field, i.e., from gluon self energy diagrams in which 
these fiavors propagate, it is expected that these logarithms will cancel against the logarithms 
present in the decay rate r(6 ^ sff) with f = u,d, s. To cancel these unphysical terms, we 
will include the Og contribution to this process in section [V|. 



V. VIRTUAL CORRECTIONS TO THE DECAY WIDTH FOR b 



sg 



We are now ready to write down the renormalized version of the matrix M''™(6 -^ sg) 
element for b ^ sg, where the virtual order a^ corrections are included. We obtain: 



M-°(6 ^sg) = ^ V^Vt, \cf + ^ 



72 



Arc 

{sg\Os{^l)\b)tree■ 



C?(£i In ^ + n) + C°(4 In ^ + r,) + 



(69) 



, m, m) and m, 



The quantities ii, ri, £2, ^2, and /g'^'^ are given in eqs. 

respectively. As eq. (^) shows, Cf^ is the only Wilson coefficient needed to NLL precision. 

For the following, it is useful to decompose it as 



a 



eff 



^8 +^(^8 



(70) 



The symbol (s(7|08(a*) |^)tree in eq. (^) denotes the tree level matrix element of Os{fi), which 
contains the running 6-quarks mass and the strong running coupling constant at the scale /x. 
In order to get expressions where the b- quark mass enters as the pole mass, and the strong 
coupling constant enters as gsirrih), we rewrite {sg\Os{lj)\b)tvee as 



{sg\Os{p)\b)tree = {sg\08\b)t 



2as . nib 4:as as ^ , mb 
1 + In -— + — /5oln — 

7T jJ, 6 IT AtT jJi 



Po 



23 



(71) 



where we made use of eqs. (|]) and ( [A18| ). The symbol {sg\Os\b)tT:ee then stands for the tree 
level matrix element of Og in which rnb{jj) and gs have to to be identified with the pole mass 
nib and gs{nib), respectively. (See also the discussion after eq. ( |i3| ) and eq. (PD). Inserting 
eqs. ( |70D and ([7l|) into eq. (^) we obtain: 



M^'^'ib -^ sg) 






Cg' +(8 + /3o) In — Cg' 



C?(4 In ^ + n) + C°(£2 In ^ + r^) + C"''^'^ fr 



/^ 



/^ 



16 

y 



•-^8 



+ 



{Sg\0s\b)tree . (72) 
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To obtain the decay width F^^''* from M'"'^°(6 -^ sg) is straightforward. We get: 



Tvirt 



2 Cts ^O,off 



247r4 



\Grv:MU{cr) +^cl 



2C8^'^f^ + 2(8 + /3o)ln^C°'^^ 



/i 



_^C°'^«^ + 2 C?(4 In ^ + Re(ri)) + 2 C^l^^ In ^ + Re(r2)) 



yO,eS 



+2C8^'^"Re(/r)(l-6) 






^2e 



l + 2e- -(7r2-16)e^ 



(73) 



We note that due to the infrared poles present in /g*^" the phase space integrations have 
been done consistently in c? = 4 — 2e dimensions, which leads to the last two extra factors 
in the last term in eq. (173|) . The other factor, (1 — e), in the last term in eq. ([73D, is due 
to the fact that all the {d — 2) possible transverse polarizations of the emitted gluon were 
taken into account. 



VI. Og CONTRIBUTION TO THE DECAY WIDTH r(6 -^ sff) 



As discussed at the end of section |V|, we should take into account the contribution of 
the operator Os to the process b -^ sff (/ = u,d,s), in order to cancel the unphysical 
logarithms of the form ln(m///i) in the virtual corrections to 6 — > sg. The Os contribution 
to the decay width T^ib — > sff) yields 



Tsib ^ sff) 



ml\GpV,:Vt,C's^''\' 



727r5 



■ a. 



In 



2m/ 



(74) 



Comparing this result with F"^* in eq. (^), we see explicitly, that the mentioned logarithms 
indeed cancel. 



VII. MATRIX ELEMENTS FOR GLUON BREMSSTRAHLUNG 



In this section we discuss the gluon bremsstrahlung corrections to 6 — *> sg, i.e. the matrix 
element for the process b -^ sgg, associated with the operators Oi, O2 and Og. For literature 
on the analogous corrections to 6 ^ 57, we refer to [^ . 



A. Bremsstrahlung associated with 0\ and O2 

We first discuss the matrix element of O2. There are two diagrams contributing; they 
are displayed in d) and e) of fig. |^. The sum of diagram d) and the one with the two gluons 
interchanged is denoted by Jap- Its analytic form is obtained by putting r^ = and r^ = 
in the expression for Jap in eq. (p]): 
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Ja^ = np{'i.r){T\T^}+f~,{q,r) 



rpA rpB 



(75) 



This expression is understood to be contracted with the polarization vectors e°'{q) and e^{r) 
of the gluons. The diagram e) in fig. 0, denoted by 5*^^, is color antisymmetric and can be 
written as 



qAB _ q- 

^0/3 ~ '-'a/3 



rpA rpB 



(76) 



where S^n reads {t = {2qr)/m 



S. 



Q/3 



_9l_ 

32772 



4 / /i \ 2^ 1 4 



3 \mcJ e 3 



8Gi(t) + 8G2(t) 



[t/9: 



aP — <iyal3 



2-fpr^ + 27„g^] L . (77) 



The functions Gj(t) {i = —1, 0, 1, ...) are defined as 



Gi{t) = dxx' ln[l - tx{l - x) - iS] . 
Jo 



(78) 



The Ward identities r^T^^ = q'^T^a = 0, stated in [Q, imply that 



T, 



a/S 



327r2 



E{a, f3, r) - E{a, (3, q) - E{I3, ^. ?) 7^ + E{a, r, q) -^ 



LAt23 . (79) 



General considerations (or a straightforward calculation which makes use of the explicit 
expressions for the functions Gi and Aij) imply the Ward identities 



r^if^p + S-p) = ; g"(f -^ + S'^) = 



(80) 



which can be used to cast (T^^ + 5^^) into the simple form 

^2 



Taf3 + S^p - 22^2 



^' ii-i) V-^-9.Alm,,. 



qr 



(81) 



To summarize, the matrix element M|""''™'^ = {sgg\02\b) can be written as 



^brems 



f+^{T^,T^} + {f-^ + S, 



a/B 



A rpB 



rpji J. 



(82) 



where T^ and (T^^ + S^^^) are given in eqs. (^) and (|81|), respectively. The functions 
Az23 and Aiu occurring in these expressions, can be written in terms of Go{t) and G_i(t) 
(t = {2qr)/m1) defined in eq. ([75|): 



Az23 = -2 -^^ ; Azi7 

6 



2 t + 6G'_i(t)-12G'o(t) 

3 t 



(83) 
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FIG. 7. Bremsstrahlung diagrams associated to Og and O2. Circle-crosses denote possible 
gluon emissions. Note that picture a) actually represents four Feynman diagrams (obtained by in- 
terchanging the gluons) and the one in d) represents two diagrams (again: including the interchange 
of the gluons). 

The explicit form of G-i(t) and Go(t) is given in appendix 0. Note that these results 
are ultraviolet finite. As the subsequent phase space integrals do not generate infrared 
singularities, it is consistent to retain terms up to order e° only in eq. (^21). 



Due to the specific color structure of the operator Oi, the diagram e) in fig. |^ does 
not contribute and the color antisymmetric part encoded in T~p is also absent. The matrix 
element M}'''''™^ = {sgg\Oi\b) is therefore proportional to T^, reading 

A, B and a, h are the color indices of the gluons and the quarks, respectively. 



B. Bremsstrahlung associated with Og 

The Feynman diagrams contributing to the matrix element M^'''^™'^ = {sgg\Oii\h) are 
shown in a), b) and c) in fig. 0. Similar to M^'''^™'^ in eq. (0), one can decompose Mg'^'^™*^ 
into a color symmetric- and a color antisymmetric part: 

The diagrams shown in b) and c) only contribute to R~a, while the diagrams in a) contribute 
to both R~^ and R^^. As the calculation of these tree level diagrams is straightforward, we 
do not give the explicit expressions for R^^ and i?~^. 



VIII. DECAY WIDTH FOR b -^ sgg 

The total matrix element M^^^'^^lb -^ sgg) can be written as 



rbrems ^ T/*1 

V2 



^Mems _ __^ y*^y^^ 
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where the three terms on the r.h.s., given in eqs. (p4D, (|82D and ([85|), correspond to the 
contributions of the operators Oi, O2 and Og, respectively. The coefficients Ci and C2 are 
understood to be the following linear combinations of the Wilson coefficients Ci and C2 
appearing in the effective Hamiltonian (j^): 

Ci = 2 C*i ; C2 = C2 — qCi ; (87) 

We note that in eq. (|86|) only the leading order pieces of the Wilson coefficients are needed. 
The expression for the decay width reads in d dimensions: 

^pbrems^^ _^ ^^^^ = J_ f (27rf S'^ip -p -q- r)TM^^^^ dfx(p')dfi(q)dfx(r) , 

ZUlh J 



where p, p', q, r are the four-momenta of the 6-quark, s-quark, and the gluons. |M'^'''=™'^|| is 
obtained by squaring the matrix element M^^'^'^^, followed by summing/ averaging over spins 
and color of the final/initial state particles. The factor (1/2) due to the two gluons in the 
final state is also absorbed there. 

The phase space integrals are plagued with infrared and collinear singularities. Configu- 
rations with one gluon flying collinear to the s-quark are regulated by a small strange quark 
mass rUs, while configurations with two collinear gluons, or one soft gluon are dimensionally 
regularized. As in the calculations of the virtual corrections, we write the dimension as 
d = 4 — 2e. (Note that e has to be negative in order to regulate the phase space integrals). 

When squaring M^'''^™'^ in eq. (|5BD, nine terms are generated, which we denote for obvious 
reasons by (61, O^;), (61, O2*), (Oi,0^), (02,0^), (Os,^^*), (02,0*), (08,0^, (08,02*), and 
(08, Og). It turns out that all terms except (08,0g) are free of infrared and collinear 
singularities. We therefore can put rris = in these terms and evaluate the phase space 
integrals in d = 4 dimensions. Denoting this finite contribution to the decay width by 
ptoems^ we get: 

The superscripts (+) and (— ) on the various r-quantities refer to color even and color odd 
contributions, respectively. The result is represented as a two dimensional integral over the 
energies Eq and Er of the gluons in the rest frame of the 6-quark. Eq and Er vary in the 
range 



Eqe 






Er e 



rrib nib 

2 '"2 



(90) 



The various r-quantities, in which all the scalar products are understood to be expressed in 



terms of Eg and Er, read: 



r+ = Cf 24 \Ai2s\^ 2ml[ml - 2(gr)] 
4, = Cn \^^23\' 2ml[ml - 2iqr)] 
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C^ 12 |Aii7p 2 [16{pqf - 16{pq){qr) - 8ml{pq) + 6m^(gr) + 



m. 



r+ = 2 Ci ^2 8 I A«23|' 2m2[m2 - 2(gr)] 
r+ = 2(71 C"'^'^ 8 Re(Az23) IGmUqr) 
4 = 2(72 C°'^«^ f Re(A^23) 16m2(gr) 



0,eflF 



T-^ = 2C2 C^'^ 12 Re(Aii7 



-Ami 



''TT'biPl) + ^(pr) — 2ml{pq)^ — 2m\{pry 



-2m1{pq){pr) + 4(pg)^(pr) + 4(pr)^(pg) /[(m)!^"")] 



(91) 



were the functions /\in and A223 are given in eq. ([53|). As these function are rather 
complicated, the integrals over Eq and E,. are done numerically. 

We now turn to the (Ogj O^) contribution, denoted by Fgg*'™^. Without going too much 
into the details, we would like to mention that some care has to be taken when summing 
over the {d — 2) transverse polarizations of the gluons. These sums are of the form 



d-2 



Y^e^^{k)er{k) = -g^'' + k^r + k''r 



(92) 



r=l 



where the vector / satisfies the condition {fk) = 1, with k being the four-momentum of 



the gluon. It turns out that both terms involving / on the r.h.s in eq. (|9^) contribute to 

)ar 
obtain (with p = {ms/rribY) 



the color antisymmetric part of Fgg'^™''. After a lengthy, but straightforward calculation, we 
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for the color symmetric part, and 
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for the color antisymmetric part. V is defined as 



V 



as ml 

247r4 
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The total decay with for b —>■ sgg is then given by 

pbrems^^ -. sgg) = T^f^^^"^' + T] 






where the three terms on the r.h.s. are given in eqs. (|89|), (|93D and 



(93) 



(94) 



(95) 



(96) 
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IX. COMBINED NLL BRANCHING RATIO FOR b ^ sg AND b -^ sgg 



In this section we combine the decay widths for the virtually corrected process b ^ sg 
and the bremsstrahlung process b —>■ sgg to the decay width, which we call r^'"'"(6 — > sg). 
We also absorb in this quantity the Og induced contribution to the process b — > sf f, where 
f = u,d, s, as discussed at the end of section |^ and in section ^ The expression for r^'''*, 



which contains the lowest order contribution to the decay width for b -^ sg, together with 
its virtual corrections, may be found in eq. (|73D . The result for the bremsstrahlung process, 
-pbrems jg gjygj^ jj^ gq (^6|). From the explicit formulas for T^"^ and r'^'''^™'^ one can see that 
the infrared singularities and those coUinear singularities, which are regulated by ejR cancel 
in the sum. The same also happens with the coUinear singularities which are regularized by 
the parameter p = {nis/mhY. The terms containing logarithms of the light quark masses 
rri/, present in the result for r^^''*, are cancelled when combined with r8(6 -^ sf f) in eq. 
(0). Putting together the individual pieces, we obtain 
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(97) 



where Fg^^™*^, given in eq. (|89|), contains all the bremsstrahlung corrections except those 
originating from the (08,0|) interference. The quantities £i, ri, £2 and r2 stem from the 
virtual corrections; they are given in eqs. (|55D, (|56D, (|5^ and (|^), respectively. On the other 
hand, fg and rg contain information from the real part of the virtual corrections, encoded 
in Re(/8'^"); the contributions from the (08,0|) interference of the gluon bremsstrahlung 
process; and the Og contribution to the process b — > sfj: The explicit expressions for £§ 
and rs (which is real by definition) read 
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We would like to stress that all scale dependent quantities in eq. (|97|) are understood to be 
evaluated at the scale /x, unless indicated explicitly in the notation. 

To prepare the discussion on the numerical size of the NLL QCD corrections, it is useful 
to cast the final result ( P?]) into another form: 



pNLL/T X _ (^s{fnb)mf^ 2 I r)|2 , pbrems 

i [b-^sg)- — ^^ — \GFVtsVtb\ H^\ +ifin 
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with 
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The modulus square of D is understood to be taken in the same way as the in the virtual 
contributions, i.e., by systematically discarding the 0{al) term. In this sense, the quantity 
D can be viewed as an effective matrix element. 

We would like to mention that £i, £2 and (£§ + 8 + [3q) are identical to the anomalous 
dimension matrix elements 7;^8' , 728^ , and 738^ , respectively. This is of course what has 
to happen: Only in this case the leading scale dependence of Cg'*' (/i) gets compensated by 
the second term in eq. ( |100| ). 

The NNL branching ratio B^^^ip -^ sg) is then obtained as 



^NLL 



B^^^{h -^ sg) 



{h -. sg) 






lexp 
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where E^^^ denotes the experimental semileptonic branching ratio of the 5-meson. T^i stands 
for the theoretical expression of the semileptonic decay width of the i?-meson. Neglecting 
non-perturbative corrections of the order (AqcD/'^b)^, Tsi reads (with Xc = {rric/mb)) 



Tib 



ceup 



1927r3 



\Vcb\^ g{xc) 



Ztc 



(102) 



where the phase space function g{xc) reads 

g{xc) = 1 - 8x1 - 24x1 \nxc + Sxl - xl ■ 



(103) 



The analytic expression for hsi{xc) can be found in ref. |2^. The approximation 

h,i{xc) = -3.341 + 4.05 (xc - 0.3) - 4.3 (x, - 0.3)^ (104) 



holds to an accuracy of 1 permille in the relevant range 0.2 < Xc < 0.4. 

We note that in the numerical analysis of B^^^{b -^ sg) we systematically expand the 
expression for the branching ratio (|101|) in a^, dropping terms of O(a^). 



A short remark concerning the LL branching ratio is in order: For the decay width 



T {b -^ sg), we use the expression 



T^^ib -. sg) 



as{fi)ml 



2An^ 



GfV;M' (C8^^'^^(/x) 



(105) 



The LL branching ratio for b ^ sg is then obtained as in eq. ( |101| ), but by discarding the 
radiative corrections in Tgi. 



X. NUMERICAL RESULTS FOR THE COMBINED BRANCHING RATIO 

Before we present the numerical result for the branching ratio B^^^{b —>■ sg), we discuss 
the sizes of the various NLL corrections at the level of the function D, defined in eq. (|100|) 
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FIG. 8. Scale (^) dependence of the function D (see eq. ( |10C1| )) in various approximations: The 
long-dashed line shows Cg'^ ; the short-dashed line corresponds to putting ri = r2 = rg = 0; the 
dotted line is obtained by only putting r2 = 0; the solid line shows the full function D. See text. 



(anticipating that the finite bremsstrahlung corrections in eq. (p9|) are relatively small). 
We already mentioned that the terms containing the explicit logarithms of the ratio {rrib/fi) 
get compensated by the scale dependence of the first term on the r.h.s. of eq. ( [1001 ). 
This feature can be observed in fig. p, when comparing the two dashed lines. The long- 
dashed line represents only the first term Cg of the function D, while the short-dashed 
line shows D, in which ri, r2 and rg are put to zero. As expected, the short-dashed line 
has a milder /i-dependence. When switching on also ri and rg (but keeping r2 = 0), the 
resulting curve, shown by the dotted line, stays close to the short-dashed curve and the 
scale dependence remains mild. However, when switching on also r2, the situation changes 
drastically. The resulting solid line, which represents the full NLL D function, implies that 
the term containing the two- loop quantity r2, induces a large NLL correction. As this large 
correction term contains a factor as{fi) C2(yu), it is of no surprise, that the NLL prediction 
for the function D suffers from a relatively large scale dependence, as illustrated by the solid 
line. 

The NLL branching ratio B^^^{b -^ sg) is then obtained as described in section |X|. The 
result is shown by the solid line in fig. ^. For the input values, we take: m^ = (4.8 ± 0.2) 
GeV, {mjmb) = 0.29 ± 0.02, a,(m^) = 0.119 ± 0.003, iV^^WKfeP = 0.95 ± 0.03, B^^^ = 
(10.49 ±0.46)%, and m^"^ = (175 ±5) GeV. As the scale dependence is rather large, we did 
not take into account the error due to the uncertainties in the input parameters. Based on 
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FIG. 9. Branching ratio B(b — > sg) as a function of the scale /i in various approximations: The 
dashed and the sohd hnes show the LL and the NLL predictions, respectively; the dotted line is 
obtained by putting ri = r2 = rg = r|j°gj^g = in the NLL expression for r^^^(6 — > sg) in eq. 
See text. 



(106) 



fig. ^ we obtain the NLL branching ratio 

B^^^ib -^ sg) = (5.0 ± LO) x lO'l 



We would like to stress that the NLL corrections drastically enhance the LL value (see 
dashed line in fig. ^ for which one obtains 

B^^{b -^ sg) = (2.2 ± 0.8) x 10"^ . 



(107) 



As already mentioned in the discussion of the function D, the main enhancement is due to 
the virtual- and bremsstrahlung corrections to b -^ sg, calculated in this paper. At the level 
of the branching ratio, this fact is illustrated by the dotted line in fig. y, which is obtained 
by discarding Fg^"™^ and by switching off ri, r2 and rg in the expression for F^'"'"(6 — ^ sg) 
(see eq. (^). 

The largest uncertainty due to the physical input parameters on B^^^{b —>■ sg) results 
from the charm quark mass. The dependence of B^^^{b — > sg) on rric is illustrated in fig. 



T0| , where Xc = rric/mb is varied between 0.27 and 0.31. Choosing ^ = nib, the resulting 



uncertainty amounts to ~ ±6%. 

XI. NUMERICAL EVALUATION OF THE CHARMLESS DECAY RATE 

In this section we investigate the impact of the NLL QCD corrections to b —>■ sg on the 
inclusive hadronic charmless decay rate of the B meson. At the quark level, we take into 
account the hadronic processes 
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FIG. 10. NLL Branching ratio B {h -^ sg) as a function of the scale ji for the three value of 
the ratio Xc = mc/rrH). See text. 



qqq 



sg-. 



(108) 



where q = d,s and q' = u, d, s. As we do not distinguish between AS* = and AS" = 
1 contributions, we can safely neglect the CKM suppressed decay mode b -^ dg. More 
precisely, we calculate the CP-averaged branching ratio 



B, 



r(&->x^) + r(6^x^) 
2rsi 
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where X^ stands for the final states listed in eq. ( |108|) . In the numerical results for B^ we 
will insert Fgi as given in eq. ( [L02| ), i.e., we do not make an a^ expansion of l/T^i in eq. 
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|109| ). The charmless hadronic decay rate B^ then reads 



B, = B. 



sg 



E 



^q'q'q 



While the 0{as) corrections to semileptonic processes have been known for a long time (see 
e.g. ref. |2^), the NLL corrections to the hadronic processes in eq. ( |108| ) with 3 quarks 
in the final state had a long history and were completed to a large extent only recently by 
Lenz et al. 0,0; however, current-current type corrections to the penguin operators are still 
missing. To briefly summarize the history, it is useful to decompose the NLL expressions for 
the decay widths of these processes into various pieces. Taking as an example the process 
b -^ uud, we write as in ref. IB! : 



r(6 — > uud) 



a. 



r(o) + :ii [AT, 



+ ALpeng + ALg + ALw] 



1111 
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The first two terms in the square bracket in eq. ( |111|) describe the effect of current 



current and penguin diagrams involving the operators^ Oi and O2. AFg hkewise contains 
the matrix element of the operator Og. The remaining part AFw of the NLL contribution 
is made of the corrections to the Wilson coefficients multiplying the tree-level amplitudes 
in T^^' . In this approximation, the matrix elements of the penguin operators Os,...,Oq only 
enter at tree level. As the expressions for the r.h.s. of eq. (|111|) are explicitly given in ref. 



1^, we do not give them here. For later reference, we denote this approximation (in lack of 
a better word) by "approxl". 

Later, in ref. 0, the same authors added the contributions of the penguin diagrams 
associated with the penguin operators to the decay matrix elements and took into account 
the interference with the tree level matrix element of the operator O2 in the decay width. In 
addition, they took into account the square of the matrix element of the penguin diagram 
associated with O2. Although being of next-to-NLL, this term is numerically relatively large, 
as it is multiplied with C|. These new contributions can be absorbed into the quantity 
AFnew, which is understood to be added to the terms in the bracket in eq. ( |111|) . As the 



extraction of AFnew from ref. is straightforward, we do not give the explicit expression. 
This approximation, which contains - up to the current-current type corrections to the 
penguin operators - the full NLL contribution to the hadronic three body decays, is called 
"approx2" . 

We note that the approximation where only the current-current type corrections AFcc 
were considered together with the shifts AFw in the Wilson coefficients has existed for a 
long time |^. We denote this approximation by "approxO" in the numerical discussion. 

In table |I| we present numerical results for the charmless hadronic branching ratio B^ 
in the various approximations mentioned above. The process b — » sg, encoded in Bgg in 
eq. (IllOl) is taken into account in the columns "approxO", "approxl" and"approx2" at LL 



precision. The last column includes in addition the NLL corrections to b ^ sg which were 
calculated in this paper. Table H was produced with the following input parameters: 



rUb = (4.8 ± 0.2) GeV, /x = rUb, {mjrrib) = 0.29 ± 0.04, 

as{mz) = 0.119 ± 0.003, m^^' = (175 ± 5) GeV, -B^f^ = (10.49 ± 0.46)% 

|Ks| = 0.22, \Vcb\ = 0.038, iKb/Kfel = 0.095 ±0.035, 5 = 60° ± 30°. (112) 

The central value for iVub/KftI corresponds to the (improved) Wolfenstein parameters p = 



0.20 and rj = 0.37 ||2^. The remaining entries of the CKM matrix are then obtained as 



described in detail in [^|. We note that the averaged charmless hadronic branching ratio is 



practically independent of 6, as already observed in ref. 0]. 

The numbers in column "approx2" are very similar to those in table 1 of ref. 0. The 
small discrepancy is due to the omission of the 1/ml power corrections in our work. 

Staring from the numbers in column "approxO", table || illustrates, that the various 
improvements shown in the other columns are relatively large, tending to increase B^. In 



^Note, that the authors of refs. PJ6[ use the old operator basis [22| 
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input 


approxO 


approxl 


approx2 


with NLL b^ sg 


as in (|112D 


1.32 


1.50 


1.62 


1.88 


H = mb/4. 


3.86 


3.21 


3.34 


3.62 


II = mb/2 


2.06 


2.09 


2.18 


2.43 


H = 2m}, 


0.96 


1.14 


1.28 


1.55 


\Vub/Vcb\ = 0.06 


0.94 


1.13 


1.24 


1.50 


\Vub/Vcb\ = 0.07 


1.03 


1.22 


1.33 


1.59 


\Vub/Vcb\ = 0.08 


1.14 


1.32 


1.44 


1.69 


\Vub/V,b\ = 0.09 


1.26 


1.44 


1.55 


1.81 


\Vub/V,b\ = 0.10 


1.39 


1.57 


1.69 


1.94 


\Vub/V,b\ =0.11 


1.54 


1.72 


1.83 


2.09 


\Vub/V,b\ = 0.12 


1.70 


1.87 


1.99 


2.25 


\Vub/Vcb\ = 0.13 


1.87 


2.05 


2.16 


2.42 


Xc = 0.25 


1.14 


1.32 


1.45 


1.69 


Xc = 0.27 


1.22 


1.41 


1.53 


1.78 


Xc = 0.29 


1.32 


1.50 


1.62 


1.88 


Xc = 0.31 


1.44 


1.61 


1.72 


1.99 


j;c = 0.33 


1.57 


1.74 


1.84 


2.12 



TABLE II. Table for the charmless hadronic branching ratio B^ (in %) in the various approxi- 
mations discussed in the text. Unless specified explicitly in the first column, the input parameters 



correspond to the central values in eq. (112). 
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particular, the NLL corrections to b ^ sg are of similar importance as the corrections 
calculated in ^,0. 

For iV^b/Kfel = 0.095 we obtain the charmless hadronic branching ratio 

B^ = {l.88t',il)%, (113) 

where the error corresponds to a variation of Xc = {mc/mf,) and of the renormalization 
scale fi in the ranges 0.25 < Xc < 0.33 and 0.5 < fj^/rrih < 2.0. The corresponding errors 
are added in quadrature. The experimental uncertainty in as{mz) has a smaller impact 



and the errors due to the remaining input parameters in eq. ( [L12| ) are negligible. The 
large renormalization scale dependence of this result is expected to be weakened once the 
current-current type corrections to the penguin operators are included. 

So far, we have considered the charmless hadronic branching ratio B^. To obtain the 
total charmless branching ratio B{B -^ no charm), one has to add twice the charmless 
semileptonic branching ratio B{B -^ X„£z7^), for £ = e and £ = jj, |]2^ (the contribution for 
£ = r, as well as radiative decay modes can be safely neglected): 

B{B ^ XJI7,) = (0.17 ± 0.03) % X (^^^) • (114) 

For iKb/Kfel = 0.095, we find 

B{B -^ no charm) = (2.22+° J°) % . (115) 

The experimental result for the total charmless branching ratio reads 

^exp^^ ^ no charm) = (0.2 ± 4.1) % , (116) 



obtained in ref. [^ from CLEO data [^ . 



XII. NUMERICAL PREDICTIONS IN THE PRESENCE OF ENHANCED Cf 

As discussed in the introduction, the theoretical prediction of the semileptonic branch- 
ing ratio and the charm multiplicity are compatible with the experimental findings if the 
renormalization scale is allowed to be as low as rrih/A. Both predictions are, however, at the 
lower side and therefore an enhancement of the charmless hadronic branching ratio B^ by 
new physics would lead to a better agreement. It is therefore still conceivable that B^ is 
considerably larger than in the standard model (SM). 

In the SM the initial conditions for C^^q and Cg are generated at a scale /i = O^mw) 
by the one-loop bsg vertex function. Due to the fact that the W^-boson only couples to left- 
handed quarks, only chromo magnetic operators proportional to mb (and m^) are generated. 
In extensions of the SM, however, also chromomagnetic operators where rrib (or m^) is 
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FIG. 11. Branching ratio B{b -^ sg) as a function of / = C| (m^yj/Cg ' (mw). For the 



exact definition of /, see eq. (117). The dotted (sohd) curve shows the LL (NLL) approximation. 
The dashed curve is obtained by switching off the matrix elements of the operators Oi and O2 ■ 



replaced by the mass of a heavy particle propagating in the loop, can be generated [^ . 
Such operators potentially lead to large contributions to 6 — *> sg. In the following we will 
perform a model independent analysis of the impact of enhanced Cg on B^, emphasizing 
the role of the NLL corrections to b -^ sg. We assume that only chromomagnetic operators 
with the same helicity structure as Os in the SM are generated which can then be described 
as a shift in Cg. For simplicity, we further assume that the CKM structure of the new 
contribution is the same as in the SM, hence neglecting the possibility of new CP-violating 
phases, by assuming the shift in Cg to be real. 



In fig 

a,/(47r 
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\n\ we investigate the impact of 
on the branching ratio for b — 






enhanced C^ (mw) = C-'g'""(mvi/) + 
sg. In the NLL approximation for this 



branching ratio, both, Cg'^ {mw) and Cg'^ (rnw) enter. In general, it is expected that the 
two pieces get different new physics shifts. For purpose of illustration, we assume however 
that both pieces are the same multiple / of the SM counterparts, i.e., we assume that 
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The dotted curve shows the LL prediction of B{b -^ sg) as a function of /, while the solid 
curve shows the NLL prediction. It is expected that for large enhancement factors, the 
matrix elements of the operators Oi and O2 become unimportant; this feature is illustrated 
by the dashed line, which is obtained by switching off these matrix elements. The NLL 
corrections (for large enhancement factors) amount to almost 50% of the LL prediction. 

In fig. |T2|, the impact of enhanced Cg on the charmless hadronic branching ratio B^ is 
illustrated. The dotted curve includes the NLL corrections to the decay modes with three 
quark in the final state and the LL result for B{b —>■ sg) (see "approx2" in section |X| ), while 
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FIG. 12. Charmless hadronic branching ratio Bj as a function of / = Cg {my\r)/C^ ' {mw 



For the exact definition of /, see eq. (117). The dotted (solid) curve includes the LL (NLL) 
approximation for B{b — > g). The NLL corrections to the decay modes with three quark in the 



final state (see "approx2" in section XI) are included in both cases. 



the solid curve also includes the NLL corrections to B{b -^ sg). For a given value of B^ 
(from an ideal measurement), Cs{rnw) can be measured in principle. To illustrate this, we 
take the hypothetical value B^ = 5%. The two solutions for the enhancement factor / are 
/ = 7 and / = —9 when using the dotted curve; including NLL corrections to b —>■ sg (solid 
curve), enhancement factors with smaller absolute values do the job, viz. / = 5 and / = —8. 



XIII. SUMMARY 

In this paper we presented a detailed calculation of the 0{as) virtual corrections to the 
decay width T{b — ^ sg). The most difficult part, the two-loop diagrams associated with the 
operators Oi and O2 which from the numerical point of view play a crucial role, was obtained 
by using Mellin-Barnes techniques. Also complete expressions for the corresponding 0{as) 
bremsstrahlung corrections to b ^ sg were given. The combined result is free of infrared 
and collinear singularities, in accordance with the KLN theorem. 

The renormalized virtually corrected matrix element {sg\Os\b) contains logarithms of the 
form \n{mf/fi) (/ = u, d, s, c, b), which for the light flavors {u, d, s) represent a special kind 
of singularity. Keeping in mind that these terms originate from the renormalization factor 
Z3 of the gluon field, i.e., from gluon self energy diagrams in which these flavors propagate, 
we argued that these singularities cancel against the logarithms present in the decay rate 
r(6 -^ sf f) with f = u,d, s. We therefore included the Og contribution to T{b -^ sf f) for 
/ = u,d,s. 
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Taking into account the existing next-to-leading logarithmic (NLL) result for the Wilson 
coefficient C^^, a complete NLL result for the branching ratio B^^^{b -^ sg) was obtained. 
Numerically, we found B^^^ = (5.0 ± 1.0) x 10~^, which is more than a factor of two larger 
than the leading logarithmic result B^^ = (2.2 ± 0.8) x 10~^. 

We then investigated the impact of these corrections on the inclusive charmless hadronic 
branching ratio B^ of S-mesons. We found that the NLL corrections calculated in this paper 
are of similar importance as NLL corrections to 6-quark decay modes with three quarks in 
the final state, which were presented by Lenz et al. [|,§. 

Finally, the impact of the NLL corrections to b —>■ sg on B^ was studied in scenarios, 
where the Wilson coefficient Cg is enhanced by new physics. For a given value of B^ (from 
an ideal measurement), Csirnw) can be measured in principle. To illustrate this, we took 
the hypothetical value B^ = 5%. The two solutions for the enhancement factor / are / = 7 
and / = —9, using the LL approximation for B{b -^ sg); including NLL corrections to 
b -^ sg, somewhat smaller enhancement factors (/ = 5 and / = —8) are needed to obtain 
the hypothetical value B^ = 5%. 
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APPENDIX A: NEXT-TO-LEADING ORDER WILSON COEFFICIENTS 



In this appendix we present the exphcit formulas which allow to calculate the Wilson 
coefficients needed in this paper. 



In section |A1| , we give the results for the Wilson coefficients at the matching scale fiw, 
which is usually taken to be of order mw Section |A2| is devoted to the the Wilson coefficients 
at the scale fib, where /x^ is of order nih. We give an explicit expression for Cf^^fih) at NLL, 
which is new. To make this appendix self-contained, we also repeat the results for the Wilson 
coefficients Ci(/i{,) and C2{fib) which are needed only to LL precision in our application. 



1. NLL Wilson coefficients at the matching scale fiw 



To give the results for the effective Wilson coefficients C^^^ at the matching scale fj,w in a 
compact form, we writeQ: 



— Oj [flw) ■ 



C^i {^^w) — Ci ' ifJ'W) + 



(Al) 



The LL Wilson coefficients at this scale are well known 33,34 . 



^0,eflf 



C2' ifJ'w) 



cr"(M = o 



C7' ifJ'w) 



X 

24 

X 



(2 = 1,3,4,5,6) 
-8x^ + 3x2 + i2x - 7 + (18x2 - 12x) Inx' 
(x-l)4 _ 

—x^ + 6x^ — 3a; — 2 — 6x In a; 
(x-lY 



(A2) 



■yO,eflf/ 



■yO.eff/ 



The coefficients Cy'^ (fiw) and Cg'*^ (fJ-w) are functions of x = mf/rri^^. Note that there 
is no explicit dependence of the matching scale fiw in these functions. Whether there is 
an implicit /xiy-dependence via the t-quark mass depends on the precise definition of this 
mass which has to be specified when going beyond leading logarithms. If one chooses to 
work with mt{fi]y), then there is such an implicit /xvy^dependence of the lowest order Wilson 
coefficient; in contrast, when working with the pole mass rrit there is no such dependence. 
We choose to express our NLL results in terms of the pole mass rrit. 

The NLL pieces Q ''^ (fJ'w) of the Wilson coefficients have an explicit dependence on the 
matching scale /ivK and for i = 7,8 they also explicitly depend on the actual definition of 
the t-quark mass. Initially, when the heavy particles are integrated out, it is convenient to 
work out the matching conditions Q '"^ (fiw) for i = 7,8 in terms of rnt{fiw)- Using eq. (H), 



^Note that Cf{fi) = Ci{p) by definition for i = 1, ...,6. 
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it is then straightforward to get the corresponding result expressed in terms of the pole mass 
mt- One obtains for i = 1, ..., 6: 



2 ' 



Cl'^'^ifxw) = ^o + ^ln 



2 i„ ^^w 



1 ' 



C'''''(M = 0, (z = 2,3,5,6) (A3) 



with 



xx^ + llx-lS a;2 4a;2 - 16x + 15 , 2, 2 ,,,, 

Eo = ^ -^ — -A ^ ^Inx Inx . (A4) 

" 12(a;-l)3 6(a;-l)4 3 3 ^ ^ 

For i = 7, 8, we split Q ''^ (/ivr) into three terms: 

C^^^(/iH.) = W-. + M, In 4^ + T, fin ^ - ^ V (A5) 



The first two terms Wi and Mj would be the full result when working in terms of mi(/ivv)- 
Ti results when expressing mt{fiy^r) in terms of the pole mass rrit in the corresponding lowest 
order coefficients. Thus, for i = 7, 8, the term Tj is obtained as 

T, = 8x ^^^'2^'^^ - (A6) 

The explicit form of the functions Wi, Mi and Tj reads 

-16a;4 - 122x3 + 80^^ - 8a; ^ / 1 \ 6x^ + 46x3 - 28^^ , „ 

Wj = z Li2 1 H z TT \n^x 

^ 9(x-l)4 'V xj 3(x-l)5 

-102x5 - 588x^ - 2262x3 + 3244x2 - 1364x + 208 , 

H ; 77 Inx 

81(x-l)5 

1646x^ + 12205x3 - 10740x2 + 2509x - 436 
^ 486(x-l)4 

-4x^ + 40x3 + 41x2 + X , / 1\ -17x3-31x2 „ 

"^ = 6(I3T)3 ^'^ ['-x) + 2(x - 1)^ '" " 

-210x5 + 1086x4 + 4893x3 + 2857x2 - 1994x + 280 , 

H ; — mx 

216(x-l)5 

737x4 _ i4io2a;3 - 28209x2 + 610x - 508 
+ 



1296(x- 1)4 

_ 82x5 + 301x4+703x3-2197x2 + 1319x-208 -(162x^ + 1242x3-756x2) Inx 
^'" 81(x-l)5 
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Ms 



T. = ^ 



T« = 2x 



77x5 _ 475^,4 _ 11113,3 ^ gQ73,2 ^ 10423; _ 140 +(918x^ + 1674x2) Inx 

108(x-l)5 

47x3 - 63x2 + 9a, + 7 _ (133,3 _^ 303,2 _ 24^,) 1^3, 



(x-l)5 

-x^ — 9x2 _^ 93, _^ 1 _^ (53.2 _^ Q^^ Ynx 
(x-l)5 



(A7) 



The dilogarithm Li2(x) is defined by 



^ rft 



Li2(x) = -/ -ln(l-t) 



(A8) 



2. NLL Wilson coefficients at the low scale fib 



The evolution from the matching scale fiw down to the low-energy scale fib is described by 
the renormalization group equation 



f^^Cf{fi) = Cf{fi)^^^{fi). 



(A9) 



The initial conditions C^^^{fiw) for this equation are given in section \KA\ , while the anomalous 
dimension matrix ■y^J^ up to order a^ can be found in ref. ^^. For completeness we display 
the result here. The anomalous dimension matrix can be expanded perturbatively as 



7/(/i)-"^^^^-°'^'' 



Air 



■7 



J* 



+ 



(47r 



1 2 %■« 



l,eff 



+ Oial 



(AlO) 



where matrix 7^^° is given by 



{^n 



/-4| 

12 







V 






2 
9 








208 
243 





4 
3 








416 
81 





52 
3 





2 


176 
81 


40 

9 


100 
9 


4 
9 


5 

6 


152 
243 





256 
3 





20 


6272 
81 


256 
9 


56 
9 


40 
9 


2 
3 


4624 
243 














32 
3 














32 
9 



173 \ 
162 ' 

70 
27 
14 

27 

587 
162 

6596 

27 

4772 

81 



28 
3 



(All) 



l,eff . 



and in the MS scheme with fully anticommuting 75, ■jjl is 
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Hn 



\ 



355 
9 

_35 
3 











502 

27 

_28 

3 











1412 
243 


1369 

243 


134 
243 


35 
162 


818 
243 


3779 
324 


416 
81 


1280 
81 


56 

81 


35 

27 


508 
81 


1841 
108 


4468 


31469 


400 


3373 


22348 


10178 


81 


81 


81 


108 


243 


81 


8158 


59399 


269 


12899 


17584 


172471 


243 


243 


486 


648 


243 


648 


251680 


128648 


23836 


6106 


1183696 


2901296 


81 


81 


81 


27 


729 


243 


58640 


26348 


14324 


2551 


2480344 


3296257 



243 








243 



243 



162 


















2187 
4688 
27 
2192 



81 



729 



4063 
27 



(A12) 



The solution of eq. ( |A9| ) , obtained through the procedure described in |2^, yields for the 
coefficient Cg^{fj,b), which we decompose as 



4:71 



Cg^^^i/^.) = cr^(/i,) + 



(A13) 



the LL term 



^0,eflf/ \ 14 Q g a'^0,off/ \ 



(A14) 



i=l 



and the NLL contribution 

Cl'^'^if^b) = r/iQ^'^^(/XH^) + 6.7441 (r/i - r^^^) C^^ifiw) 



(A15) 



i=l 



The symbol r] is defined as ?7 = as{\i-w)las{iih)\ the vectors a,, a'j ,/i^, e'j, //, k\ and l\ read 
{a,} = 1 1|^ 1|^ A^ _1|^ 0.4086, -0.4230, -0.8994, 0.1456} 
{a^} = { if, 0.4086, -0.4230, -0.8994, 0.1456} 
{/i^} = { §§§, -0.9135, 0.0873, -0.0571, 0.0209} 



{e'J = {2.1399, 0, 0, 0, -2.6788, 0.2318, 0.3741, -0.0670} 
{/;} = {-5.8157, 0, 1.4062, -3.9895, 3.2850, 3.6851, -0.1424, 0.6492} 
{A;^} = {3.7264, 0, 0, 0, -3.2247, 0.3359, 0.3812, -0.2968} 
{Q = {0.2169, 0, 0, 0, -0.1793, -0.0730, 0.0240, 0.0113} . (A16) 

As already mentioned earlier, we neglect the contributions of the operators 03,...,0q in 
our analysis for b —^ sg, as their Wilson coefficients are rather small. We therefore only list 
the results for the coefficients C]^^(/ib) and C2^{fii,), which are needed to LL precision only: 
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(A17) 



When calculating NLL results in the numerical analysis, we use the NLL expression for 
the strong coupling constant: 



as(/i) 






1- 



(3i asimz) lnt;(/i) 
/3o 47r f(/i) 



(A18) 



with 



v{fi) = 1 - /?o ^^ In I 



(A19) 



where /?o = 3 and /5i = -3^ (for 5 flavors). However, for LL results we always use the LL 
expression for as(/i), i.e., /3i is put to zero in eq. ( |A18| ). 



APPENDIX B: ONE-LOOP FUNCTIONS G_i(t) AND Go{t) 



In this appendix we give the explicit results for the functions G'_i(t) and Golt) needed 
in eq. (PB]). Evaluating the integral in eq. (^) for z = — 1, 0, one obtains: 



G-i{t) 



_zf + 2 In^ (^^±1^) - 2z7r In (^^±|^) ; 
-4 - 2 arctan^ (y^) + 27r arctan (y^ 



t>4 
< t <4 



(Bl) 



Go(t) 



-2 + 2w^ln( ^+f^ 



?7ri 






—2 — 2a/ ^arctan 



t>4 



(\/¥)+-\/¥; 0<t<4 



(B2) 
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